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CLAIMS 



What is claimed is: 



1. 



A method for automated system identification comprising: 



selecting a model structure; 

generating at least one reference signal for input into a system; 

retrieving a plurality of input signals and a plurality of output signals 
from said system; 

performing system identification on said model structure using said 
plurality of input signals, said plurality of output signals, and said at least one 
reference signal to obtain a point model; and 

verifying accuracy of said point model. 

2. The method according to claim 1, further comprising: 
calculating a cost vector for said model structure; and 

selecting a model order based on said cost vector associated with said 
model structure. 

3. The method according to claim 2, wherein said model structure 
includes at least one model parameter. 

4. The method according to claim 1, further comprising: 
providing at least one operating condition for said system; 
providing a sampling frequency and a frequency bandwidth covered by 

said model structure; and 
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defining a plurality of identification experiments according to said at 
least one operating condition, said sampling frequency, and said frequency 
bandwidth. 

5. The method according to claim 4, wherein said model structure 
includes a plurality of experimental parameters determined by said plurality of 
identification experiments. 

6. The method according to claim 1, further comprising verifying 
linearity of said system and detecting non-linear manifestations of said system. 

7. The method according to claim 1, further comprising storing said 
at least one reference signal into a reference table and storing said plurality of 
input signals and said plurality of output signals into an input/output table. 

8. The method according to claim 1, wherein said model structure is 
selected from a group consisting of finite impulse response (FIR), 
autoregressive with external input (ARX), autoregressive moving average with 
external input (ARMAX), autoregressive moving average (ARMA), 
autoregressive autoregressive with external input (ARARX), autoregressive 
autoregressive moving average with external input (ARARMAX), output error 
(OE), Box-Jenkins (BJ), and Ordinary Differential Equations (ODE). 

9. The method according to claim 4, further comprising: 
performing each identification experiment of said plurality of 

identification experiments in said system; and 
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obtaining said plurality of input signals and said plurality of output 
signals from said each identification experiment. 

10. The method according to claim 1, wherein said reference signal is 
selected from a group consisting of a chirp signal, a pseudo random binary 
sequence, a sum of sinusoids, and a wavelet. 

11. The method according to claim 1, wherein said reference signal is 
generated to obtain a large output signal-to-noise ratio and to guarantee a 
linear operation regime for said point model. 

12. The method according to claim 1, further comprising: 
automatically detecting at least one outlier value in said plurality of 

output signals; and 

removing said at least one outlier value from said plurality of output 
signals. 

13. The method according to claim 12, further comprising replacing 
said at least one outlier value with a predetermined value calculated using a 
filter. 

14. The method according to claim 12, wherein said detecting further 
comprises: 

building a filter using said plurality of input signals and said plurality of 
output signals; 

computing said at least one outlier value using said filter; 
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comparing said at least one outlier value with a predetermined threshold 
value; and 

storing said at least one outlier value if said at least one outlier value is 
greater than said predetermined threshold value. 

15. The method according to claim 14, wherein said detecting 
requires a plurality of iterations, each iteration being related to a time value. 

16. The method according to claim 1, wherein said performing further 
comprises identifying an input/ output model and a disturbance model within 
said point model. 

17. The method according to claim 16, wherein said input/output 
model is unstable and said disturbance model is determined using said 
input/output model. 

18. The method according to claim 16, further comprising calculating 
an input/output uncertainty parameter within said input/output model, and 
calculating a disturbance uncertainty parameter within said disturbance model. 

19. The method according to claim 1, wherein said verifying further 
comprises analyzing whether a plurality of innovation signals, derived from 
said plurality of output signals, are white stochastic signals uncorrected with 
past measurements. 

20. The method according to claim 2, wherein said verifying further 
comprises analyzing said plurality of input signals and said plurality of output 
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3 signals retrieved using a value of zero for said at least one reference signal and 

4 increasing said model order to account for unrepresented system dynamics. 

1 21. The method according to claim 18, wherein said verifying further 

2 comprises analyzing frequency regions corresponding to said input/ output 

3 uncertainty parameter and said disturbance uncertainty parameter and 

4 modifying said at least one reference signal by increasing a sweep time 

5 corresponding to said frequency regions. 

1 22. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for automated system identification comprising: 

4 selecting a model structure; 

5 generating at least one reference signal for input into a system; 

6 retrieving a plurality of input signals and a plurality of output signals 

7 from said system; 

8 performing system identification on said model structure using said 

9 plurality of input signals, said plurality of output signals, and said at least one 

10 reference signal to obtain a point model; and 

1 1 verifying accuracy of said point model. 

1 23. The computer readable medium according to claim 22, wherein 

2 said method further comprises: 

3 calculating a cost vector for said model structure; and 

4 selecting a model order based on said cost vector associated with said 

5 model structure. 
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1 24. The computer readable medium according to claim 23, wherein 

2 said model structure includes at least one model parameter. 

1 25. The computer readable medium according to claim 22, wherein 

2 said method further comprises: 

3 providing at least one operating condition for said system; 

4 providing a sampling frequency and a frequency bandwidth covered by 

5 said model structure; and 

6 defining a plurality of identification experiments according to said at 

7 least one operating condition, said sampling frequency, and said frequency 

8 bandwidth. 

1 26. The computer readable medium according to claim 25, wherein 

2 said model structure includes a plurality of experimental parameters 

3 determined by said plurality of identification experiments. 

1 27. The computer readable medium according to claim 22, wherein 

2 said method further comprises verifying linearity of said system and detecting 

3 non-linear manifestations of said system. 

1 28. The computer readable medium according to claim 22, wherein 

2 said method further comprises storing said at least one reference signal into a 

3 reference table and storing said plurality of input signals and said plurality of 

4 output signals into an input /output table. 

1 29. The computer readable medium according to claim 22, wherein 

2 said model structure is selected from a group consisting of finite impulse 
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3 response (FIR), autoregressive with external input (ARX), autoregressive 

4 moving average with external input (ARMAX), autoregressive moving average 

5 (ARMA), autoregressive autoregressive with external input (ARARX), 

6 autoregressive autoregressive moving average with external input 

7 (ARARMAX), output error (OE), Box-Jenkins (BJ), and Ordinary Differential 

8 Equations (ODE). 

1 30. The computer readable medium according to claim 25, wherein 

2 said method further comprises: 

3 performing each identification experiment of said plurality of 

4 identification experiments in said system; and 

*H 5 obtaining said plurality of input signals and said plurality of output 

6 signals from said each identification experiment. 
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1 31. The computer readable medium according to claim 22, wherein 

2 said reference signal is selected from a group consisting of a chirp signal, a 

3 pseudo random binary sequence, a sum of sinusoids, and a wavelet. 

1 32. The computer readable medium according to claim 22, wherein 

2 said reference signal is generated to obtain maximum output signal-to-noise 

3 ratio and to guarantee a linear operation regime for said point model. 

1 33. The computer readable medium according to claim 22, wherein 

2 said method further comprises: 

3 automatically detecting at least one outlier value in said plurality of 

4 output signals; and 
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removing said at least one outlier value from said plurality of output 
signals. 

34. The computer readable medium according to claim 33, wherein 
said method further comprises replacing said at least one outlier value with a 
predetermined value calculated using a filter. 

35. The computer readable medium according to claim 33, wherein 
said detecting further includes: 

building a filter using said plurality of input signals and said plurality of 
output signals; 

computing said at least one outlier value using said filter; 

comparing said at least one outlier value with a predetermined threshold 
value; and 

storing said at least one outlier value if said at least one outlier value is 
greater than said predetermined threshold value. 

36. The computer readable medium according to claim 35, wherein 
said detecting requires a plurality of iterations, each iteration being related to a 
time value. 

37. The computer readable medium according to claim 22, wherein 
said performing further comprises identifying an input/output model and a 
disturbance model within said point model. 
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1 38. The computer readable medium according to claim 37, wherein 

2 said input/output model is unstable and said disturbance model is determined 

3 using said input/output model. 

1 39. The computer readable medium according to claim 37, wherein 



2 said method further comprises calculating an input/output uncertainty 

3 parameter within said input/output model, and calculating a disturbance 

4 uncertainty parameter within said disturbance model. 

1 40. The computer readable medium according to claim 22, wherein 

2 said verifying further comprises analyzing whether a plurality of innovation 

a "3 

4s 3 signals, derived from said plurality of output signals, are white stochastic 

y ?. 

flj 4 signals uncorrected with past measurements. 



m 1 41. The computer readable medium according to claim 23, wherein 

JU 2 said verifying further comprises analyzing said plurality of input signals and 

r; 3 said plurality of output signals retrieved using a value of zero for said at least 

~ 4 one reference signal and increasing said model order to account for 

P 5 unrepresented system dynamics. 

1 42. The computer readable medium according to claim 39, wherein 

2 said verifying further comprises analyzing frequency regions corresponding to 

3 said input/ output uncertainty parameter and said disturbance uncertainty 

4 parameter and modifying said at least one reference signal by increasing a 

5 sweep time corresponding to said frequency regions. 
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1 43. An article of manufacture comprising a program storage medium 

2 readable by a computer and tangibly embodying at least one program of 

3 instructions executable by said computer to perform a method for automated 

4 system identification, said method comprising: 

5 selecting a model structure; 

6 generating at least one reference signal for input into a system; 

7 retrieving a plurality of input signals and a plurality of output signals 

8 from said system; 

9 performing system identification on said model structure using said 

10 plurality of input signals, said plurality of output signals, and said at least one 

1 1 reference signal to obtain a point model; and 

12 verifying accuracy of said point model. 

1 44. The article of manufacture according to claim 43, wherein said 

2 method further comprises: 

3 calculating a cost vector for said model structure; and 

4 selecting a model order based on said cost vector associated with said 

5 model structure. 

1 45. The article of manufacture according to claim 44, wherein said 

2 model structure includes at least one model parameter. 

1 46. The article of manufacture according to claim 43, wherein said 

2 method further comprises: 

3 providing at least one operating condition for said system; 

4 providing a sampling frequency and a frequency bandwidth covered by 

5 said model structure; and 
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6 defining a plurality of identification experiments according to said at 

7 least one operating condition, said sampling frequency, and said frequency 

8 bandwidth. 

1 47. The article of manufacture according to claim 46, wherein said 

2 model structure includes a plurality of experimental parameters determined by 

3 said plurality of identification experiments. 

1 48. The article of manufacture according to claim 43, wherein said 

2 method further comprises verifying linearity of said system and detecting non- 

3 linear manifestations of said system. 

1 49. The article of manufacture according to claim 43, wherein said 

2 method further comprises storing said at least one reference signal into a 

3 reference table and storing said plurality of input signals and said plurality of 

4 output signals into an input/ output table. 



W 1 50. The article of manufacture according to claim 43, wherein said 

2 model structure is selected from a group consisting of finite impulse response 

3 (FIR), autoregressive with external input (ARX), autoregressive moving average 

4 with external input (ARMAX), autoregressive moving average (ARMA), 

5 autoregressive autoregressive with external input (ARARX), autoregressive 

6 autoregressive moving average with external input (ARARMAX), output error 

7 (OE), Box-Jenkins (BJ), and Ordinary Differential Equations (ODE). 

1 51. The article of manufacture according to claim 46, wherein said 

2 method further comprises: 
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performing each identification experiment of said plurality of 
identification experiments in said system; and 

obtaining said plurality of input signals and said plurality of output 
signals from said each identification experiment. 

52. The article of manufacture according to claim 43, wherein said 
reference signal is selected from a group consisting of a chirp signal, a pseudo 
random binary sequence, a sum of sinusoids, and a wavelet. 

53. The article of manufacture according to claim 43, wherein said 
reference signal is generated to obtain maximum output signal-to-noise ratio 
and to guarantee a linear operation regime for said point model. 

54. The article of manufacture according to claim 43, wherein said 
method further comprises: 

automatically detecting at least one outlier value in said plurality of 
output signals; and 

removing said at least one outlier value from said plurality of output 
signals. 

55. The article of manufacture according to claim 54, wherein said 
method further comprises replacing said at least one outlier value with a 
predetermined value calculated using a filter. 

56. The article of manufacture according to claim 54, wherein said 
detecting further includes: 



34 




3 building a filter using said plurality of input signals and said plurality of 

4 output signals; 

5 computing said at least one outlier value using said filter; 

6 comparing said at least one outlier value with a predetermined threshold 

7 value; and 

8 storing said at least one outlier value if said at least one outlier value is 

9 greater than said predetermined threshold value. 

1 57. The article of manufacture according to claim 56, wherein said 

2 detecting requires a plurality of iterations, each iteration being related to a time 

3 value. 

O 

tSsr 

nj 1 58. The article of manufacture according to claim 43, wherein said 

s n i: 

-p 2 performing further comprises identifying an input/ output model and a 

m 

Ul 3 disturbance model within said point model. 



?* 1 59. The article of manufacture according to claim 58, wherein said 

jf 2 input/ output model is unstable and said disturbance model is determined 

S3 3 using said input/output model. 

1 60. The article of manufacture according to claim 58, wherein said 

2 method further comprises calculating an input/ output uncertainty parameter 

3 within said input/ output model, and calculating a disturbance uncertainty 

4 parameter within said disturbance model. 

1 61. The article of manufacture according to claim 43, wherein said 

2 verifying further comprises analyzing whether a plurality of innovation signals, 



35 



derived from said plurality of output signals, are white stochastic signals 
uncorrelated with past measurements. 

62. The article of manufacture according to claim 44, wherein said 
verifying further comprises analyzing said plurality of input signals and said 
plurality of output signals retrieved using a value of zero for said at least one 
reference signal and increasing said model order to account for unrepresented 
system dynamics. 

63. The article of manufacture according to claim 60, wherein said 
verifying further comprises analyzing frequency regions corresponding to said 
input/ output uncertainty parameter and said disturbance uncertainty 
parameter and modifying said at least one reference signal by increasing a 
sweep time corresponding to said frequency regions. 

64. A method for automated system identification comprising: 
qualifying a system; 

performing an identification experiment procedure on said system to 
obtain a plurality of input signal values and a plurality of output signal values; 
filtering said plurality of output signal values to obtain point model data; 

and 

validating a point model obtained using said point model data. 

65. The method according to claim 64, wherein said qualifying further 
comprises: 

calculating a cost vector associated with each model structure of a 
plurality of model structures for said system; 
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selecting one model structure based on said associated cost vector; and 
selecting a model order based on said one model structure and said 
associated cost vector. 

66. The method according to claim 65, wherein said calculating 
further comprises calculating said cost vector as a function of a risk of local 
minima factor, characteristic to said system, a computational cost factor, 
characteristic to said model structure, and an equipment time factor related to a 
number of identification experiments necessary to obtain said point model. 

67. The method according to claim 65, further comprising 
transmitting said cost vector to a real-time planner module for selection of said 
one model structure. 

68. The method according to claim 65, further comprising 
transmitting said cost vector to a user for selection of said one model structure. 

69. The method according to claim 65, further comprising 
transmitting said cost vector to a processing module for selection of said one 
model structure. 

70. The method according to claim 65, wherein said one model 
structure selected further includes at least one model parameter. 

71. The method according to claim 65, wherein said qualifying further 
comprises: 

providing at least one operating condition for said system; 
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4 providing a sampling frequency and a frequency bandwidth covered by 

5 said one model structure selected; and 

6 defining a plurality of identification experiments according to said at 

7 least one operating condition, said sampling frequency, and said frequency 

8 bandwidth. 

1 72. The method according to claim 71, wherein said one model 

2 structure selected further includes at least one experimental parameter 

3 determined by said plurality of identification experiments. 

1 73. The method according to claim 64, wherein said qualifying further 

41 2 comprises: 

Ul 

3 performing at least one qualification test on said system; and 

4 calculating a qualification vector based on said at least one qualification 

5 test. 

1 74. The method according to claim 73, wherein said performing said 

2 at least one qualification test further comprises: 

3 injecting a first reference signal value into said system to obtain a first 

4 output signal value; 

5 injecting a second reference signal value, obtained by scaling said first 

6 reference signal value by a predetermined factor, to obtain a second output 

7 signal value; and 

8 comparing said first output signal value to said second output signal 

9 value to verify linearity of said system and to detect non-linear manifestations 
10 of said system. 
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1 75. The method according to claim 65, wherein said model structure 

2 is selected from a group consisting of finite impulse response (FIR), 

3 autoregressive with external input (ARX), autoregressive moving average with 

4 external input (ARMAX), autoregressive moving average (ARMA), 

5 autoregressive autoregressive with external input (ARARX), autoregressive 

6 autoregressive moving average with external input (ARARMAX), output error 

7 (OE), Box- Jenkins (BJ), and Ordinary Differential Equations (ODE). 
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1 76. The method according to claim 73, further comprising deciding 

2 whether said system is qualified based on results from said at least one 

3 qualification test and terminating said qualifying if said results are outside of a 

4 predetermined range. 
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1 77. The method according to claim 71, wherein said performing 

2 further comprises: 

3 generating at least one reference signal value for input into said system; 

4 performing each identification experiment of said plurality of 

5 identification experiments in said system using said at least one reference signal 

6 value; and 

7 obtaining said plurality of input signal values and said plurality of 

8 output signal values from said each identification experiment performed. 

1 78. The method according to claim 77, wherein said performing 

2 further comprises storing said plurality of input signal values into an input 

3 storage device and said plurality of output signal values into an output storage 

4 device. 



1 79. The method according to claim 77, wherein said generating said at 

& * 2 least one reference signal value is based on at least one model parameter 

3 associated with said model structure and at least one experimental parameter 

4 associated with said model structure and determined by said plurality of 

5 identification experiments. 

1 80. The method according to claim 77, wherein said generating said at 

2 least one reference signal value produces maximum output signal-to-noise ratio 

3 and guarantees a linear operation regime for said point model. 

1 81. The method according to claim 77, wherein said at least one 

2 reference signal includes at least one signal selected from a group consisting of 
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3 a chirp signal, a pseudo random binary sequence, a sum of sinusoids, and a 

4 wavelet. 

1 82. The method according to claim 77, wherein said performing 

2 further comprises storing said at least one reference signal value into a 

3 reference storage device. 

1 83. The method according to claim 78, wherein said generating 

2 further comprises: 

3 retrieving one output signal value from said output storage device; 

4 dividing said one output signal value by said at least one reference signal 

5 value to obtain an input/ output gain; and 

6 dividing a predetermined output signal level by the input/output gain 

7 to obtain a new reference signal value. 



* 1 84. The method according to claim 77, wherein said generating is 

O 

W 2 iterative, being performed repetitively for each identification experiment of said 

n *** 

□ 3 plurality of identification experiments. 

O 

1 85. The method according to claim 64, wherein said filtering further 

2 comprises: 

3 automatically detecting at least one outlier value in said plurality of 

4 output signal values; and 

5 removing said at least one outlier value from said plurality of output 

6 signal values. 
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86. The method according to claim 85, wherein said detecting further 
comprises: 

constructing a filter using said plurality of input signal values and said 
plurality of output signal values; 

computing said at least one outlier value as a difference between a 
predetermined output signal value corresponding to said filter and one output 
signal value of said plurality of output signal values; 

comparing said at least one outlier value with a predetermined threshold 
error value; and 

storing said at least one outlier value if said at least one outlier value is 
greater than said predetermined threshold value. 

87. The method according to claim 86, wherein said filtering further 
comprises replacing said at least one outlier value with said predetermined 
output signal value calculated using said filter. 

88. The method according to claim 85, wherein said detecting is 
iterative, being performed repetitively for a plurality of time values if said at 
least one outlier value is lower than a predetermined threshold value. 

89. The method according to claim 86, wherein said comparing is 
automatically performed by a real-time planner. 

90. The method according to claim 86, wherein said comparing is 
performed by a user. 
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91. The method according to claim 86, wherein said comparing is 
automatically performed by a processing module. 

92. The method according to claim 85, further comprising: 
removing said plurality of output signal values if said at least one outlier 

value being removed exceeds a predetermined percentage of said plurality of 
output signal values; and 

iteratively performing said identification experiment procedure and said 
filtering to obtain new point model data. 

93. The method according to claim 64, wherein said validating further 
comprises: 

analyzing whether a plurality of innovation signal values, derived from 
said plurality of output signal values, correspond to a plurality of white 
stochastic signal values; and 

iteratively performing said identification experiment procedure and said 
filtering to obtain new point model data if said plurality of innovation signal 
values do not correspond to said plurality of white stochastic signal values. 

94. The method according to claim 64, wherein said validating further 
comprises: 

generating at least one reference signal value for input into said system; 

analyzing said plurality of input signal values and said plurality of 
output signal values retrieved using a zero value for said at least one reference 
signal value; 

calculating output spectral estimates for said plurality of output signal 
values; 
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9 calculating input spectral estimates for said plurality of input signal 

10 values; 

1 1 calculating a transfer function estimate as a ratio of said output spectral 

12 estimates and said input spectral estimates; 

13 comparing said transfer function estimate with said point model data; 

14 and 

15 iteratively performing said identification experiment procedure and said 

16 filtering to obtain new point model data if features of said output spectral 

17 estimates and said input spectral estimates are not present in said point model 

18 data. 

1 95. The method according to claim 64, further comprising performing 

2 an identification on said point model. 

1 96. The method according to claim 95, wherein said performing of 

2 said identification further comprises: 

3 identifying an input/output model within said point model, said 

4 input/output model being characterized by an input/ output transfer function; 

5 identifying a disturbance model within said point model, said 

6 disturbance model being characterized by a disturbance transfer function; 

7 assessing stability of said input/output model; and 

8 calculating said disturbance transfer function based on said stability of 

9 said input/output model. 

1 97. The method according to claim 96, wherein said input/output 

2 model is unstable and said calculating further comprises: 
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3 processing said input/ output transfer function to obtain at least two 

4 stable transfer functions; 

5 calculating a prediction error associated with said system based on said 

6 at least two stable transfer functions; and 

7 calculating said disturbance transfer function using said prediction error 

8 and said at least two stable transfer functions for a predetermined model 

9 structure. 

1 98. The method according to claim 97, wherein said predetermined 

2 model structure is selected from a group consisting of finite impulse response 

3 (FIR), autoregressive with external input (ARX), autoregressive moving average 

4 with external input (ARMAX), autoregressive moving average (ARMA), 
J 5 autoregressive autoregressive with external input (ARARX), autoregressive 

6 autoregressive moving average with external input (ARARMAX), output error 

7 (OE), Box-Jenkins (BJ), and Ordinary Differential Equations (ODE). 



1 99. A method for qualification of a system comprising: 

2 calculating a cost vector associated with each model structure of a 

3 plurality of model structures for said system; 

4 selecting one model structure based on said associated cost vector; and 

5 selecting a model order based on said one model structure and said 

6 associated cost vector. 

1 100. The method according to claim 99, wherein said calculating 

2 further comprises: 

3 calculating said cost vector as a function of a risk of local minima factor, 

4 characteristic to said system; 
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calculating a computational cost factor, characteristic to said model 
structure; and 

calculating an equipment time factor related to a number of 
identification experiments necessary to obtain said point model. 

101. The method according to claim 99, further comprising 
transmitting said cost vector to a real-time planner module for selection of said 
one model structure. 

102. The method according to claim 99, further comprising 
transmitting said cost vector to a user for selection of said one model structure. 

103. The method according to claim 99, further comprising 
transmitting said cost vector to a processing module for selection of said one 
model structure. 

104. The method according to claim 99, wherein said one model 
structure selected further includes at least one model parameter. 

105. The method according to claim 99, further comprising: 
providing at least one operating condition for said system; 
providing a sampling frequency and a frequency bandwidth covered by 

said one model structure selected; and 

defining a plurality of identification experiments according to said at 
least one operating condition, said sampling frequency, and said frequency 
bandwidth. 
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1 106. The method according to claim 105, wherein said one model 

2 structure selected further includes at least one experimental parameter 

3 determined by said plurality of identification experiments. 

1 107. The method according to claim 99, further comprising: 

2 performing at least one qualification test on said system; and 

3 calculating a qualification vector based on said at least one qualification 

4 test. 

1 108. The method according to claim 107, wherein said performing said 

2 at least one qualification test further comprises: 

3 injecting a first reference signal value into said system to obtain a first 



yi 4 output signal value; 
\7* 5 injecting a second reference signal value, obtained by scaling said first 



6 reference signal value by a predetermined factor, to obtain a second output 

7 signal value; and 

p a 8 comparing said first output signal value to said second output signal 

^ 9 value to verify linearity of said system and to detect non-linear manifestations 

0 10 of said system. 



1 109. The method according to claim 99, wherein said model structure 

2 is selected from a group consisting of finite impulse response (FIR), 

3 autoregressive with external input (ARX), autoregressive moving average with 

4 external input (ARMAX), autoregressive moving average (ARMA), 

5 autoregressive autoregressive with external input (ARARX), autoregressive 

6 autoregressive moving average with external input (ARARMAX), output error 

7 (OE), Box-Jenkins (BJ), and Ordinary Differential Equations (ODE). 



47 



110. The method according to claim 107, further comprising deciding 
whether said system is qualified based on results from said at least one 
qualification test and terminating said qualifying if said results are outside of a 
predetermined range. 
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111. A method for performing an identification experiment on a 
system comprising: 

generating at least one reference signal value for input into said system; 

performing each identification experiment of a plurality of identification 
experiments in said system using said at least one reference signal value; and 

obtaining a plurality of input signal values and a plurality of output 
signal values from said each identification experiment performed. 

112. The method according to claim 111, further comprising storing 
said plurality of input signal values into an input storage device and said 
plurality of output signal values into an output storage device. 

113. The method according to claim 111, wherein said generating said 
at least one reference signal value is based on at least one model parameter 
associated with a model structure of a plurality of model structures associated 
with said system and at least one experimental parameter associated with said 
model structure and determined by said plurality of identification experiments. 

114. The method according to claim 111, wherein said generating said 
at least one reference signal value produces maximum output signal-to-noise 
ratio and guarantees a linear operation regime for said system. 

115. The method according to claim 111, wherein said at least one 
reference signal includes at least one signal selected from a group consisting of 
a chirp signal, a pseudo random binary sequence, a sum of sinusoids, and a 
wavelet. 
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116. The method according to claim 111, further comprising storing 
said at least one reference signal value into a reference storage device. 

117. The method according to claim 112, wherein said generating 
further comprises: 

retrieving one output signal value from said output storage device; 

dividing said one output signal value by said at least one reference signal 
value to obtain an input/output gain; and 

dividing a predetermined output signal level by the input/ output gain 
to obtain a new reference signal value. 

118. The method according to claim 111, wherein said generating is 
iterative, being performed repetitively for each identification experiment of said 
plurality of identification experiments. 

119. A method for filtering a plurality of output signal values obtained 
for a system comprising: 

automatically detecting at least one outlier value in said plurality of 
output signal values; and 

removing said at least one outlier value from said plurality of output 
signal values. 

120. The method according to claim 119, wherein said detecting 
further comprises: 

constructing a filter using a plurality of input signal values and said 
plurality of output signal values; 
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computing said at least one outlier value as a difference between a 
predetermined output signal value corresponding to said filter and one output 
signal value of said plurality of output signal values; 

comparing said at least one outlier value with a predetermined threshold 
error value; and 

storing said at least one outlier value if said at least one outlier value is 
greater than said predetermined threshold value. 

121. The method according to claim 120, further comprising replacing 
said at least one outlier value with said predetermined output signal value 
calculated using said filter. 

122. The method according to claim 119, wherein said detecting is 
iterative, being performed repetitively for a plurality of time values if said at 
least one outlier value is lower than a predetermined threshold value. 

123. The method according to claim 120, wherein said comparing is 
automatically performed by a real-time planner. 

124. The method according to claim 120, wherein said comparing is 
performed by a user. 

125. The method according to claim 120, wherein said comparing is 
automatically performed by a processing module. 

126. The method according to claim 119, further comprising: 
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removing said plurality of output signal values if said at least one outlier 
value being removed exceeds a predetermined percentage of said plurality of 
output signal values; and 

providing a second plurality of output signal values for further 
processing. 

127. A method for validating a point model obtained for a system 
comprising: 

generating at least one reference signal value for input into said system; 

performing at least one identification experiment in said system using 
said at least one reference signal value; 

obtaining a plurality of input signal values and a plurality of output 
signal values from said at least one identification experiment performed; 

analyzing a plurality of innovation signal values, derived from said 
plurality of output signal values; and 

validating accuracy of said point model if said plurality of innovation 
signal values corresponds to a plurality of white stochastic signal values. 

128. The method according to claim 127, further comprising: 
analyzing said plurality of input signal values and said plurality of 

output signal values retrieved using a zero value for said at least one reference 
signal value; 

calculating output spectral estimates for said plurality of output signal 
values; 

calculating input spectral estimates for said plurality of input signal 
values; and 
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9 calculating a transfer function estimate as a ratio of said output spectral 

10 estimates and said input spectral estimates; 

1 1 comparing said transfer function estimate with said point model; and 

12 validating accuracy of said point model if features of said output spectral 

13 estimates and said input spectral estimates are present in said point model 

1 129. A method for performing identification on a point model for a 

2 system comprising: 

3 identifying an input/ output model within said point model, said 

4 input/output model being characterized by an input/output transfer function; 

5 identifying a disturbance model within said point model, said 

O 6 disturbance model being characterized by a disturbance transfer function; 
U1 7 assessing stability of said input/ output model; and 

8 calculating said disturbance transfer function based on said stability of 

9 said input/ output model. 
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1 130. The method according to claim 129, wherein said input/ output 

2 model is unstable and said calculating further comprises: 

3 processing said input/output transfer function to obtain at least two 

4 stable transfer functions; 

5 calculating a prediction error associated with said system based on said 

6 at least two stable transfer functions; and 

7 calculating said disturbance transfer function using said prediction error 

8 and said at least two stable transfer functions for a predetermined model 

9 structure. 
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1 131. The method according to claim 130, wherein said predetermined 

2 model structure is selected from a group consisting of finite impulse response 

3 (FIR), autoregressive with external input (ARX), autoregressive moving average 

4 with external input (ARMAX), autoregressive moving average (ARMA), 

5 autoregressive autoregressive with external input (ARARX), autoregressive 

6 autoregressive moving average with external input (ARARMAX), output error 

7 (OE), Box-Jenkins (BJ), and Ordinary Differential Equations (ODE). 

1 132. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for automated system identification comprising: 

4 qualifying a system; 

5 performing an identification experiment procedure on said system to 

6 obtain a plurality of input signal values and a plurality of output signal values; 

7 filtering said plurality of output signal values to obtain point model data; 

8 and 

9 validating a point model obtained using said point model data. 

1 133. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for qualification of a system comprising: 

4 calculating a cost vector associated with each model structure of a 

5 plurality of model structures for said system; 

6 selecting one model structure based on said associated cost vector; and 

7 selecting a model order based on said one model structure and said 

8 associated cost vector. 
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1 134. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for performing an identification experiment on a system comprising: 

4 generating at least one reference signal value for input into said system; 

5 performing each identification experiment of a plurality of identification 

6 experiments in said system using said at least one reference signal value; and 

7 obtaining a plurality of input signal values and a plurality of output 

8 signal values from said each identification experiment performed. 

1 135. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for filtering a plurality of output signal values obtained for a system 

4 comprising: 

5 automatically detecting at least one outlier value in said plurality of 

6 output signal values; and 

7 removing said at least one outlier value from said plurality of output 

8 signal values. 



1 136. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for validating a point model obtained for a system comprising: 

4 generating at least one reference signal value for input into said system; 

5 performing at least one identification experiment in said system using 

6 said at least one reference signal value; 

7 obtaining a plurality of input signal values and a plurality of output 

8 signal values from said at least one identification experiment performed; 
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9 analyzing a plurality of innovation signal values, derived from said 

10 plurality of output signal values; and 

1 1 validating accuracy of said point model if said plurality of innovation 

12 signal values corresponds to a plurality of white stochastic signal values. 

1 137. A computer readable medium containing executable instructions 

2 which, when executed in a processing system, cause said system to perform a 

3 method for performing identification on a point model for a system comprising: 

4 identifying an input/output model within said point model, said 

5 input/output model being characterized by an input/output transfer function; 

6 identifying a disturbance model within said point model, said 

7 disturbance model being characterized by a disturbance transfer function; 

8 assessing stability of said input/output model; and 

9 calculating said disturbance transfer function based on said stability of 
10 said input/output model. 
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